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I. INTRODUCTION 

What is the origin of mass? This is a long-standing question and still remains unsolved 

[ 1 ]- 

The recent discovery of the Higgs particle [2, 3] may hint how to go beyond the standard 
model (SM). The measured Higgs mass and top quark mass [ 1] are such that the SM remains 
perturbative below the Planck scale [5-7]. According to Bardeen [8], “the SM does not, by 
itself, have a fine-tuning problem”. Because the Higgs mass term is the only term, which 
breaks scale invariance at the Lagrangian level in the SM, we may ask about the origin of this 
mass term. Mostly scale invariance is hardly broken by quantum anomaly [9]. Therefore, a 
dimensional transmutation can occur at the quantum level, which can be used to generate 
a la Coleman-Weinberg [10] the Higgs mass term in a classically scale invariant extension of 
the SM [11]- [ ]. Dynamical chiral symmetry breaking [46, 47] can also be used [48]-[ ]. 
The idea is the same as that of technicolor model [ >5, 56], where the only difference is that 
we now allow the existence of fundamental scalars. 

In this paper we consider the latter possibility, in particular the model studied in [ 18-52] . 
In this model the scale, generated in a QCD-likc hidden sector, is transmitted to the SM 
sector via a real SM singlet scalar S to trigger spontaneous breaking of electroweak (EW) 
gauge symmetry [48, 49] (see also [ ]). Moreover, due to the dynamical chiral symmetry 

breaking in the hidden sector there exist Nambu-Goldstone (NG) bosons, which are massive, 
because the coupling y of S with the hidden sector fermions breaks explicitly chiral symme¬ 
try. Therefore, the mass scale of the NG bosons, which are dark matter (DM) candidates, 
is not independent (as it is not the case in the most of DM models); it is smaller than 
the hidden sector scale, which is in the TeV region unless the coupling y is very small, i.e. 

< o(io- 4 ). 

As in [ 1, 52] we employ the Nambu-Jona-Lasinio (NJL) theory [46, 47] as an low-energy 
effective theory of the hidden sector and base our calculations on the self-consistent mean 
field (SCMF) approximation [ 38, 59] of the NJL theory, which is briefly outlined in Sect. III. 
In [ 1, 52] the maximal global flavor symmetry SU(3)v (along with a U( l)y) has been 
assumed. In this paper we relax this assumption and consider in detail the cases, in which 
SU(3)y is broken into its subgroups. We find in Sect. IV that the consistent parameter space 
can be considerably extended if SU(3)v is broken to its subgroup SU(2)y x U( 1)^,. The 
main reason is that, if SU(3)v is broken, a new mechanism for the DM annihilation, inverse 
DM conversion, becomes operative at finite temperature: A pair of lighter DM particles 
annihilate into a pair of heavier (would-be) DM particles, which subsequently decay into 
SM particles (mainly into two ys). 

Before we discuss the DM phenomenology of the model, we develop an effective theory 
for DM interactions (a linear sigma model) in the framework of the SCMF approximation 
of the NJL theory. Using the effective theory we compute the DM relic abundance and 
analyze the direct and indirect DM detection possibilities in Sect. IV. Sect. V is devoted 
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to Conclusion, and in Appendix A we give explicitly the NJL Lagrangian in the SCMF 
approximation in the case that SU(3)v is broken into U(l)§, x U( l)g. In Appendix B the 
inverse DM (mesons for QCD) propagators and also how the NJL parameters are fixed can 
be found. The one-loop integrals that are used in our calculations are collected in Appendix 
C. 


II. THE MODEL 

We consider a classically scale invariant extension of the SM studied in [48-52] 1 which 
consists of a hidden SU{ii c )h gauge sector coupled via a real singlet scalar S to the SM. 
The hidden sector Lagrangian of the model is written as 

£ h = -|Tr F 2 + Tr + g H YG, + g'QYB, - yS) tJj , (1) 

where G^ is the gauge field for the hidden QCD, B At is the f/(l)y gauge field, i.e. 

= cos 9wAfj, — sin dwZ^ , g' — e/ cos^w , (2) 

and the rif (Dirac) fermions {i = 1,..., rif) in the hidden sector belong to the fundamental 
representation of SU(71 c )h■ The trace in (1) is taken over the flavor as well as the color 
indices. The hidden fermions carry a common U (l)y charge Q, implying that they contribute 
only to Llyy of the gauge boson self-energy diagrams so that the S, T, U parameters remain 
unchanged. The £sm+s part of the total Lagrangian Ct = £h + £sm+s contains the SM 
gauge and Yukawa interactions along with the scalar potential 

Dsm+5 = a h(H^H) 2 + - 1 -\ hs S 2 {H^H) , (3) 

where H T = (H + , (h + iG)y/2) is the SM Higgs doublet field, with H + and G as the 
would-be Nambu-Goldstone fields 2 . The basic mechanism to trigger the EW symmetry 
breaking is very simple: The non-perturbative effect of dynamical chiral symmetry breaking 
in the hidden sector generates a robust scale which is transferred into the SM sector through 
the real singlet S. Then the mass term for the Higgs potential is generated via the Higgs 
portal term in (3), where the in front of the positive A hs is an assumption. 

A. Global Symmetries 

The Yukawa coupling of the hidden fermions with the singlet S breaks explicitly chiral 
symmetry. Therefore, in the limit of the vanishing Yukawa coupling matrix the global 
symmetry SU(nf)L x SU(rif)R x U( l)y x U(1)a is present at the classical level, where 
1/(1 )a is broken by anomaly at the quantum level down to its discrete subgroup Y 2n/ , and 
the unbroken U{l)y ensures the conservation of the hidden baryon number. The non-abelian 
part of the chiral symmetry SU{rif )r x SU{rif)R is broken dynamically down to its diagonal 

1 See also [60]. 

2 This classically scale invariant model is perturbatively renormalizable, and the Green’s functions are 
infrared finite [61, 62]. 
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subgroup SU(rif)v by the non-vanishing chiral condensates ( 0 , 00 , implying the existence 
of iTj — 1 NG bosons 0 a (a — 1,..., rij — 1). In the rif — 3 case the NG bosons are like the 
mesons in the real hadron world: 

7T° = 03 , 7^ = (0i =F i0 2 )/v / 2 , 

= ( 0 4 T *0 5 )/v / 2 , /?°(W°) = (06 + (—)*0 7 )/v / 2 , t / 8 = 08 , (4) 

where 7) 8 will mix with ?) 0 to form the mass eigenstates fj and fj '. (The'should avoid the 
confusion with the real mesons 7T° etc.) 

In the presence of the Yukawa coupling the chiral symmetry is explicitly broken; this 
is the only coupling which breaks the chiral symmetry explicitly. Because of this coupling 
the NG bosons become massive. An appropriate chiral rotation of 0* can diagonalize the 
Yukawa coupling matrix: 

Vij = Vi$ij (yi > 0 ) ( 5 ) 

can be assumed without loss of generality, which implies that U{l) ni ~ l corresponding to 
the elements of the Cartan subalgebra of 577 (n/) are unbroken. We assume that none of 
Hi vanishes so that all the NG bosons are massive. If two are the same, say y\ = y 2 , 
one 17(1) is promoted to an SU( 2 ). Similarly, if three y t s are the same, a product group 
1/(1) x 1/(1) is promoted to an SU( 3), and so on. In addition to these symmetry groups, 
there exists a discrete Z 4 , 

Z 4 : 0* ->• (exp f ( 71 / 2 ) 75 ) 0j = 7750 * and 5 ->■ -S . (6) 

This discrete symmetry is anomalous for odd rif, because the chiral transformation in (6) 
is an element of the anomalous U(1)a- If u/ is even, then the chiral transformation is an 
element of the anomaly-free subgroup Y 2n/ of U(1)a- Needless to say that this Z 4 is broken 
by a non-vanishing vacuum expectation value (VEV) of S, which is essential to trigger the 
EW gauge symmetry breaking. 


B. Dark Matter Candidates 

The NG bosons, which arise due to the dynamical chiral symmetry breaking in the hidden 
sector, are good DM candidates, because they are neutral and their interactions with the 
SM part start to exist at the one-loop level so that they are weak. However, not all NG 
bosons can be DM, because their stability depends on the global symmetries that are in 
tact. In the following we consider the case for rif = 3, which can be simply extended to an 
arbitrary rif. For rif = 3 there are three possibilities of the global symmetries: 


i) U( 1)b' x U(l) s if y\ 0 y 2 0 y 3 , 

( 7 ) 

ii) SU(2) v xU{l) s if yi = y 2 0 2/3 , 

( 8 ) 

iii) SU(3) v if y x = y 2 = y 3 , 

( 9 ) 


where we have suppressed U(l)y which always exists, and the case (iii) has been treated 
in detail in [52]. Without loss of generality we can assume that the elements of the Cartan 


5 


Table I. The NG bosons and DM candidates for n/ = 3. 



7T° 

7T + 

7f 

K° 

K+ 

k- R° 

fj 8 

U (l)y charge 

0 

0 

0 

0 

0 

0 0 

0 

B' 

0 

2 

-2 

-1 

1 

-1 1 

0 

B 

0 

0 

0 

3 

3 

-3 -3 

0 

SU{ 2) v 


3 


2 

2 

1 

SU( 3)y 
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subalgebra corresponding to U(l)§, and U(l)^ are 


B' = 

( 1 0 o\ 
0-10 

, B = 

( 1 0 0 \ 
0 1 0 


o 

o 

o 


^0 0 -2 y 


( 10 ) 


In Table I we show the NG bosons for nf — 3 with their quantum numbers. As we can see 
from Table I the NG bosons n° and fj 8 are unstable for the case (i) and in fact can decay 
into two 7 s, while for the case (ii) only fj 8 is unstable. Whether the stable NG bosons can 
be realistic DM particles is a dynamical question, which we will address later on. 


C. Perturbativity and stability of the scalar potential at high energy 

Before we discuss the non-perturbative effects, we consider briefly the perturbative part 
at high energies, i.e. above the scale of the dynamical chiral symmetry breaking in the 
hidden sector. As explained in the Introduction, it is essential for our scenario of explaining 
the origin of the EW scale to work that the scaler potential is unbounded below and the 
theory remains perturbative (no Landau pole) below the Planck scale. So, we require: 

4vr > X H , A s >0, 4tt > X HS > 0 , \y\ 2 < 4tt, (11) 

2 a / X h X s — X H s > 0 . ( 12 ) 

In the following discussion we assume that the perturbative regime (of the hidden sector) 
starts around qo = 1 TeV and gjj(qo)/An ~ 1. Although in this model the Higgs mass 
depends mainly on two parameters, A h and A hs, lowering Xu{qo) < 0.13 will destabilize the 
Higgs potential while increasing Xn(qo) > 0.14 will require a larger mixing with S, which is 
strongly constrained. Therefore, we consider the RG running of the couplings with Xn(qo) 
fixed at 0.135 and rely on one-loop approximations. In the case that the hypercharge Q of 
the hidden fermions is different from zero, these fermions contribute to the renormalization 
group (RG) running of the U( l)y gauge coupling considerably. We found that Q ^ 0.8 
should be satisfied for g' to remain perturbative below the Planck scale. 

Because of (12) the range of A s is constrained for a given A hs and A h- The larger Xhs 
is, the larger A s has to be. But there is an upper limit for A s(Qo) because of perturbativity 
. In Fig. 1 we show the allowed area in the X$ — y plane for different values of Xhs{%) with 
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Figure 1. Stability constraint. The allowed area in the A s — y plane for different values of A hs(qo) with 
Xh{qo) fixed at 0.135 (q 0 = 1 TeV) is shown, where we have used Q = 1/3 and assumed the SU(3)y flavor 
symmetry defined in (9). The green circles, red circles and blue points stand for Xhs{qo) = 0.1,0.06 and 
0 . 02 . 

Xh(qo) fixed at 0.135 in the SU(3)v case (9), i.e. y = y i = 2/2 = 2/3 3 - The green circles, 
red circles and blue points stand for Xhs(Qo) = 0.1, 0.06 and 0.02. There will be no allowed 
region for A hs(Qo) ~ 0.12. We have used Q = 1/3, but the allowed area does not depend 
very much on Q as long as Q < 0.8 is satisfied (which ensures perturbativity of the U(l)y 
gauge coupling). If SU(3)v is broken, then the vertical axis in Fig. 1 represents the largest 
among yiS. 


III. NAMBU-JONA-LASINIO METHOD 


A. NJL Lagrangian in a mean-field approximation 


Following [ ] we replace the high energy Lagrangian Ch in (1) by the NJL Lagrangian 
£njl = Tr + g'Q^B^ - yS)ip + 2 G Tr + G D (det $ + h.c.) , (13) 


where 

n^ — 1 

Ty = A(i - 75 )Vh = \ A i* Tr ^ A “( 1 - 75)^ , 

a=0 


(14) 


and A a (a = 1,... ,rij — 1) are the Gell-Mann matrices with A 0 = \J 2/3 1. The effective 
Lagrangian £njl has three dimensional parameters G,Gd and the cutoff A, which have 
canonical dimensions of —2, —5 and 1, respectively. Since the original Lagrangian Ch has 
only one independent scale, the parameters G,Gd and A are not independent. We restrict 


3 The same analysis has been performed in [52], but without including the constraint (12). 
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ourselves to n c = nj = 3, because in this case these parameters, up-to an overall scale, 
can be approximately fixed from hadron physics [58, 59]. The six-fcrmi interaction in (13) 
is present due to chiral anomaly of the axial U(1)a and is invariant under Z e , so that the 
NJL Lagrangian (13) has the same global symmetry as the high energy Lagrangian (1). 
Furthermore, as we mentioned in Sect. II A, we can assume without loss of generality that 
the Yukawa coupling matrix y is diagonal (see (5)). To deal with the non-renormalizable 
Lagrangian (13) we employ [51] the SCMF approximation which has been intensely studied 
by Hatsuda and Kunihiro [58, 59] for hadron physics. The NJL parameters for the hidden 
QCD is then obtained by the upscaling of the actual values of G,Gd and the cutoff A from 
QCD hadron physics. That is, we assume that the dimensionless combinations 

G 1/2 A = 1.82 , (~G d ) 1/5 A = 2.29 , (15) 

which are satisfied for hadrons, remain unchanged for a higher scale of A. 


Below we briefly outline the SCMF approximation. We go via a Bogoliubov-Valatin 
transformation from the perturbative vacuum to the “BCS” vacuum, which we simply denote 
by |0). This vacuum is so defined that the mesons (mean fields) are collected in the VEV 
of the chiral bilinear: 

ip = mi - 7 5 )A>|0) = (diag(<7i, d 2 , d 3 ) + f(A a ) T 0 a ) , (16) 

where we denote the pseudo NG boson fields after spontaneous chiral symmetry breaking by 
c f > a . The dynamics of the hidden sector creates a nonvanishing chiral condensate (O^^lO) 
which is nothing but — (dj)/4G. The actual value of (d*) can be obtained through the 
minimization of the scalar potential, as we describe shortly. In the SCMF approximation 
one splits up the NJL Lagrangian (13) into the sum 

£njl = + £/ , (17) 


where Cj is normal ordered (i.e. (0|£/|0) = 0), and £ 0 contains at most fermion bilinears 
which are not normal ordered. At the non-trivial lowest order only £ 0 is relevant for the cal¬ 
culation of the effective potential, the DM mass and the DM interactions. The explicit form 
for Cq can be found in Appendix A. The effective potential can be obtained by integrating 
out the hidden fermion fields in the BCS vacuum. At the one-loop level we find 


L7 


NJL 




i=l,2,3 


(18) 


*= 1 , 2,3 


where I\/(rii) is given in Eq. (Cl), and the constituent fermion masses Mj are given by 


G 


Afj <7i + yiS g^2 <^i+l^j+2i 


(19) 


where d 4 = di and a 5 = a 2 - Once the free parameters of the model ?/j, X H , X H s , X$ are given, 
the VEVs of dj and S can be determined through the minimization of the scalar potential 
ksM+s + Vnjl, where Vsm+s is defined in (3). After the minimum of the scalar potential is 
fixed, the mass spectrum for the CP-even particles h , S and d as well as the DM candidates 
with their properties are obtained. 




B. The value of y and hidden Chiral Phase Transition 


The Yukawa coupling in (1) violates explicitly chiral symmetry and plays a similar role as 
the current quark mass in QCD. It is well known that the nature of chiral phase transition 
in QCD depends on the value of the current quark mass. Therefore, it is expected that the 
value of y strongly influences the nature of the chiral phase transition in the hidden sector, 
which has been confirmed in [ ]. The hidden chiral phase transition occurs above the EW 

phase transition, where the nature of the EW phase transition is not known yet. In the 
following discussions, we restrict ourselves to 

y & 0.006 , ( 20 ) 

because in this case the hidden chiral phase transition is a strong first order transition [51] 
and can produce gravitational wave back ground [63, 64], which could be observed by future 
experiments such as Evolved Laser Interferometer Space Antenna (eLISA) experiment [65]. 
Needless to say that the smaller is y, the better is the NJL approximation to chiral symmetry 
breaking. 


IV. DARK MATTER PHENOMENOLOGY 

A. Dark Matter Masses 

Our DM candidates are the pseudo NG bosons, which occur due to the dynamical chiral 
symmetry breaking in the hidden sector. They are CP-odd scalars, and their masses are 
generated at one-loop in the SCMF approximation as the real meson masses, where we here, 
too, restrict ourselves to n c = rif = 3. Therefore, their inverse propagators can be calculated 
in a similar way as in the QCD case, which is given in Appendix B. 

First we consider the SU(3)v case (9) to obtain the DM mass ?udm 4 and the mass of 
the singlet ms for 0.001 ^ yi = 7/2 = 2/3 ~ 0.006. In Fig. 2 (left) we show the area in the 
ms~m dm plane, in which we obtain a correct Higgs mass, while imposing the perturbativity 
(11) as well as stability (12) constraints. The upper limit of tudm for a given ms is due 
to the upper limit of the Yukawa coupling (see (20)), while its lower limit comes from the 
lower limit of the Yukawa coupling, which is taken to be 0.0005 here. The upper limit for 
ms is dictated by the upper limit of A 5 , which is fixed by the perturbativity and stability 
constraints (11) and (12). The lowest value of ms, 250 GeV, comes from the lowest value 
of A s, which is set at 0.05 here. If SU( 3)y is only slightly broken, the DM mass will not 
change very much. 

We next consider the U( 1)x U (1)^ case (7). We may assume without loss of generality 
that the hierarchy y\ < y 2 < 2/3 is satisfied. In Fig. 2 (right) we show the ratio m^o/m^± 
versus 2/1 / 2 / 2 , where we have fixed y\ and 2/3 at 0.002 and 0.006, respectively. We can 


4 Since SU{2>)y is unbroken, all the DM particles have the mass which is denoted by here. 



9 



Figure 2. Left: The area in the mj-niDM plane for 0.0005 < yi = 2/2 = 2/3 < 0.006 in the SU(3)v case 
(9), where todm = m *0 = m*± = m^ 0 = m^± = ra*j. The ratio ms /ffiDM = 2 is satisfied on the blue 
dashed line, on which the resonance condition in the s-channel diagram for the DM annihilation (Fig. 10) 
is satisfied. Right: The ratio m^o/m^± versus yi/ 1/2 in the U(l)^, x U(l)g case (7), where y\ and 2/3 are 
fixed at 0.002 and 0.006, respectively. The constraints imposed on X H , X H g and X$ are such that we obtain 
a correct Higgs mass and the perturbativity (11) as well as stability (12) constraints are satisfied. 


conclude from Fig. 2 (right) that 7 T° is the lightest among the pseudo NG bosons and the 
ratio does not practically depend on the scalar couplings Xhj^hs and A 5 . The 

SU{2) V case (8) can be realized if two of are the same. There are two independent 
possibilities: (a) iji — y 2 < V 3 , and (b) t/i < t /2 = 2 / 3 - The mass spectrum for the case (b) 
is similar to that for the x U(l)g case. In particular, 7T° is the lightest among the 

pseudo NG bosons. As for the case (a) the mass hierarchy 

mj = mj 0 = mj± < m r < = mxo = rn k± < m fi (21) 

is always satisfied. 

The different type of the DM mass spectrum will have an important consequence when 
discussing the DM relic abundance. 


B. Effective interactions for DM Decay and Annihilations 

As discussed in Sect. II B, if the SU(3)v flavor symmetry is broken to U(l)g, x U(l)g, 
there will be two real decaying would-be DM particles fj, 7 r°, and three pairs of complex DM 
particles (K°, K °), K ± and 7 r ± . Here we will derive effective interactions for these DM fields 
by integrating out the hidden fermions at the one-loop order. The one-loop integrals and 
their lowest order expressions of expansion in the external momenta in the large A limit are 
given in Appendix B. Except for the 0-0-7 and 0 - 7-7 interactions, we assume SU(2)y flavor 
symmetry, i.e. 

(cf) = ( 0 * 2 )j Mi = M-2 , Zx = Zx± = Z K0 , = Zn± = Zft 0 , (22) 

where Z s are the wave function renormalization constants given in (B8). This is because, 
we have to assume at least SU(2) V for a realistic parameter space as we will see. 
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• 0 - 0 - 7 

The corresponding one-loop diagram is shown in Fig. 3, where the right diagram in Fig. 3 
yields zero contribution. 

= A^ v (gK+K-' y9 i j,K + d u K + Gx0K0yd^K°d u K 0 + G^+^-^d^7t + d u 7t ^ 

+Z^ (G K+K - z dpK+d v K- + G KOR o z d^k°d u R 0 + GWz^tt+cU-) , (23) 

where A(Z) flu = d il A(Z) u — d v A(Z) tx . The effective couplings in the large A limit are 

Gk+k-^ = 2Z R ±n c eQ ( 1 — ^^-(cr 2 ) J /^ 2 7 (M 3 ,Mi) , (24) 

G K o R o-y = 2Z Ro n c eQ ^1 — ^ 2 7 (M 3 ,M 2 ) , 

G 7r+7r - 7 = 2 Z„±n c eQ ( 1 - ^(cr 3 ) j /^ 7 (M 2 ,Mi) , 

Gk+K-Z = —twGx+K-'y , G k o R o Z = —twG K o R 0 1 , G n + n -z = —t M /G 7r + 7r - 7 , 

^ = (sin 6 ^/cos #^/) 2 — 0.3, and /^(m^mj) is given in (C4). In the SU(2) V limit, we 
obtain G K o#o 7 = Gk+k - 7 and G 7r + 7r - 7 = 0 , because /^ 7 (m fl ,m(,) —> (mb — m a )/( 487 r 2 m^) 
as mb —> m a . 


•0 —> 77 

The diagram in Fig. 4 shows the decay of fj, 7 r° and S' into two 7 s, but they can also decay 
two Zs and 7 and Z, if the processes are kinematically allowed. Using the NJL Lagrangian 
(A 2 ) and (C5) in Appendix C we find that the effective interaction takes the form 


= \v^ 

+ Ve^ ( -G 
4 V 2 


2 Gf\ 7 2 ■^■fH'A OL ^ + GrjjZ + ^Gr/Z 2 Z(_u/Z a f3 


t 0 7 2 + G n 0ryZ Ai^Zap + —G n oz^ Zfj,uZ a /3 J , 


(25) 


where in the large A limit 

0 ■ ^ (2(<r2> ■ <<r3)) ) Mfi 
+ (l - |4(2 (,ti> - {a,))) Mp - (2 - | 0 ( 2 < ff ,} _ fo))) M 3 - 1 ] , 

G^jZ twG^l, G^Z 2 ^]yG 7 -Jry2 ? 

G»v = zU«Aq 2 (1 - w) (Mf 1 - Mp) , 

G-]fOry Z t\yG-ft 0 7 2 , G 7v O Z 2 tyyG'ftO^ . 

As we see from (27), the if 0 —> 77 decay vanishes in the SU(2) V limit, because M x = 

this limit. 


G 2 


= Z 1/2 


a 


n r 


a/37T 


■<5 2 


(26) 

(27) 


M -2 in 


The decay of S into two 7 , two Z and 7 Z can be described by 

£ s ,. = s (\g S i , + G S7 Z + k:; sz . 


( 28 ) 
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where we find from (D 8 ) 

G w = yQ 2 Y, ViMi 1 , (29) 

i=l,2,3 

Gs 1 z = —twGs'y 2 , Gsz 2 = t'wGs'y 2 • (30) 

•Dark matter conversion 

The diagrams in Figs. 5 and 6 are examples of DM conversion, in which two incoming 
DM particles are annihilated into a pair of two DM particles which are different from the 
incoming ones. There are DM conversion amplitudes, which do not vanish in the SU(3)y 
limit, and those which vanish in the limit. Except the last r]K 2 Ti interaction term, the 
effective interaction term below do not vanish the SU(3)y limit. 


1 


£*4 = -G V 2 K 2 f, 2 (A' U A' U + K+K-) + -G^2 ff ( -(ft u ) 2 + tt+tt 

+gw (k°k° + k + k -) (4 (ft 0 ) 2 + n + r 
+G vK 2 7T r) f(k°k° - k + k~)n° + V2k°k~TT + + V2k°k + r 


(31) 


where 

G r 


— „ Zf,Z i/flr ( 1 


3 

+4 1 1 


G 


D 


8 G 2 


( 01 ) 


1 - 


G 


^( 2 (a 1 )-(a 3 ))j 


G 


D 


-2 1 


8 G 2 
G d 
8 G 2 


(ai)) 

G d 




-Zr n Z knc 


G 


D 


Grj^ir^ ^Z^Z^Tlr ( 1 




8 G 2 
G d 


gG2 (2 ( ai )-(a 3 ))j ip(M lt M 3 ) 


(32) 


8 G 2 




G 


I™(Mi) 

(4/Jf (M,) - /JF(M 3 )) , 

Gd 


~Z^nJ GD 


/c 2 7r 2 


+4Z*Z^n c 


G,*% = -^zfz k ZZ‘n c I 1 - ^L( ai ) 


8 G 2 
G d 
8 G 2 


(33) 


G 


D 


8 G 2 


( 01 ) 


ijkiMyM,) , 


8 G 2 


<a 3 >) I 2 /(MyM 3 ) 


(34) 


V 2 


73 71 

x 


G 


D 


8 G 2 


8 G 2 ^ 


1 “ ^+‘>) +( M >- M 3) - (l - ^(2 (<ti> - (<r 3 )) J I?(M u M 3 ) 


-^zl l 2 Z k zl , 2 n c 


G 


D 


8 G 2 


(/]f(M„M 3 )-/|f(M I )) , 


(35) 


and /I 4 (m a ,rrib) etc. are defined in (C8)-(C12) in Appendix C. We have not included the 
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contributions from the diagram like one in Fig. 6, because they are negligibly suppressed in 
a realistic parameter space, in which SU(3)v is only weakly broken. Similarly, G v k 2 w , too, 
is negligibly small ( G r)K 2 w /G 11 2 K 2 ~ 10 -4 ), so that we will not take into account the rjK 2 ir 
interactions in computing the DM relic abundance. 

•Dark matter coupling with S 

The diagrams in Figs. 7 and 8 show dark matter interactions with the singlet S. The DM 
coupling with S (Fig. 7) can be described by 

£(/> 2 S — S + +G K 2 S (K°K° + K + K ) + G n 2 S (-(7T°) 2 + 7T + 7T )\ . (36) 


Using (C15) - (C17) in Appendix C we find in the large A limit 

G d " 2 


Gmps q Zfi'n c 


— Z~n r 
3 ' 


%i ~~ ^^(^i) i) + 2 y 3 ^1 


G 


D 


8 G 2 




g^w) rj4(A/ s 

(37) 


Gk 2 s — =-2Zf>n c ( 1 


G 


D 


8 G 2 


{Vi) ) {yill 2 s(M 3 , MO + y 3 1^2 s (Mi, Ms)) 


~‘2Zj>n c 


G 


D 


8 G 2 


yiipsiMi), 


G„'i Q — —AZ^Ilr I 1 


G 


D 


(°3> ) yi7j4( M i) - 2 Z*n c 


G 


D 


y^s( M s) 


8G 2 '~°'J ai ‘ rsv '‘ I/ - “ 7r ' c \8G 2 
The DM coupling with two Ss (Fig. 8) can be described by 


= 0 2 ++Ga'»^(A- 0 a- 0 + k*k~) + G, Is .(i(* 0 ) 2 + #+t-; 


(38) 

(39) 

(40) 


where 


Gri 2 S 2 g ZfjTl c 


G 


2 /l ( 1 — ^pK 2 ^) _ ( a3 )) ) ( 2 ^5 2 (-^i) + WO) 


G 


1 7 

A/fjTlQ 

O 


8G 2 

Gd 

8G 2 


+2yl (1 - ) {U 2 MMs) + 1^2(M 3 )) 


(4^7^ s 1 (M 1 )- ! /3 2 /^ ss (M3)) , 


(41) 


Gk 2 s 2 — —^Z^n c ( 1 


G 


D 


8 G 2 


(Vi)) (y‘i IpAM. 3 , Ml) + y 2 4^(Mi, M 3 ) + y,y 3 /]^(M 1 , M 3 )) 


'^A' n c ^1 I<i> 2 S 2 (Ml) , 


(42) 


G-^2^2 2Z^n c ( 1 


G 


D 


8 G 2 


(° 3 ) ) y 2 (M|4 2 (Mi) + 7|? S 2 (MO) 


Gi 


- z * n °• ( 43 ) 

•Dark matter coupling with two ys 

The diagram in Fig. 9 shows the annihilation of 7r ± pair into two ys, where the annihilations 
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K+J<° 


K~,K° 



K+,K° 


K-.K° 



Figure 3. The 7 coupling (charge radius), which vanishes in the SU(3)v limit. The tt-tt -7 coupling- 
vanishes in the SU{2)v limit. 


into 7 Z, two Z s and also into two Ss are also possible if they are kinematically allowed. 

= \ A ^ (V) 2 + G k ^(K°R° + K+K~) + G^ 1 -^ 0 ) 2 + tt+tT) 

+ \a^Z^ (g^z 1 -^) 2 + G K 2 iZ (K°K° + K+K-) + G^ z ( l -y) 2 + tt+tT)^) 

+ \z^Z^ (c^ifj) 2 + G K 2 Z 2(K°k° + K + K~) + G n 2 Z 2^(7r 0 ) 2 + tt+tT)) , 

where A{Z) IW = d ll A{Z) lJ — d v A(Z)^. Using the approximate form (C29) and (C30) we find 

(45) 


a 


a 


G-i p'y 2 Zjfi c Q .4^(77) — ZfjTi c Q *4^(77) ) 


7T 


7 r 


G'ljZ'yZ twG'q 2 ^y 2 , G T i 2 Z 2 t^yGrq2^2^ 


a 


a 


G k 2 7 2 = Z k n c -Q 2 41 ^( 77 ) ~ Z k n c -Q AR 77 ) , 


7T 


7T 


GkZ'jz — —UvGk 2 7 2 , G K 2 Z 2 — t\ v G K ‘ 2 1 i i 


a 


G n 2~2 Zjj-Ti c Q -4s-(77) > 


7T 


G-ff2ry Z t'W/G'ft 2^2, G-ft 2 Z 2 t]yG-ft2ry2 , 


where 


A (77) = 3 


- 1 - 


G 


D 


8 G 2 


(a 3 ) M x + 


—mt 1 

8 G 2 3 


(46) 


(47) 


and . 4 ^( 77 ) = * 4 ^( 77 ) = * 4 ^( 77 ) in the SU(3)v limit. In a realistic parameter space for 
the SU(2)v case, the ratio *4*( 77 )/A?( 77 ), for instance, is at most 1.004. 

In the following discussions we shall use the effective interaction terms derived above to 
compute the DM relic abundance as well as the cross sections for the direct and indirect 
detections of DM. 


C. Relic Abundance of Dark Matter 

The SU(3)v case (9) has been discussed in [ 1, 52], and so we below consider only 
the (i) and (ii) cases, which are defined in (7) and (8), respectively. In a one-component 
DM system, the velocity-averaged annihilation cross section (va) should be ~ 10 -9 GeV -2 
to obtain a realistic DM relic abundance Vth 2 ~ 0.12. A rough estimate of the velocity- 
averaged annihilation cross section for DM conversion (Fig. 5 ) shows (va(fjfj — > 7t + n~)) ~ 
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7 



7 


5 


7 



7 


Figure 4. Decay of DM and S into two 7 s. In the SU(2) V limit 7 r° does not decay. 



Figure 5. Examples for DM conversion 

10~ 5 (1 — mf/m-) 1 / 2 GeV -2 , where it vanishes if SU(3)v is unbroken. The reason for the 
large annihilation cross section for DM conversion is that the coupling of the hidden fermions 
to the hidden mesons is of 0(1): There is no coupling constant for the coupling as one can 
see from the NJL Lagrangian (A2). That is, the annihilation cross section for DM conversion 
is about four orders of magnitude larger than that in an ordinary case, unless the masses of 
the incoming and outgoing DMs are almost degenerate. 


1. (i)U(l) s ,xU(l) s 

There exists a problem for the U( 1)^, x U( l)g case (7), which we will discuss now. As 
we have found in the previous subsection, the lightest NG boson in the U( 1)^, x V(\) t 
case is always the lightest between the neutral ones within the one-loop analysis in the NJL 



Figure 6. This DM conversion vanishes in the SU( 3)y limit. 
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Figure 7. DM coupling with one S. In a realistic parameter space there is an accidental cancellation 
between these two diagrams. 



Figure 8. DM coupling with two S s. These diagrams contribute to the DM relic abundance if the mass 
of S is comparable with or less than the DM masses. 

approximation and that without of loss of generality we can assume it is j f°. Its dominant 
decay mode is into two ys. The decay width can be calculated from the effective Lagrangian 
(25): 

r(ft° ->• 77 ) = 9Z Q^ 3 a mlo ^1 - ( 1 /Ml - 1 /M 2 ) 2 

~ 10” 6 x Q 4 ml o (1/Mi - 1/M 2 ) 2 , (48) 

which should be compared with the expansion rate H of the Universe at T = o, 

r 7 77 ^ ~ 7 x 10 8 g 4 (m i o/M 1 ) 3 (AM/M 1 ) 2 [T eV/M^ + 0(AM 3 ) , (49) 

H 

where AM = Mi — M 2 . Therefore, unless the U(1)y charge Q of the hidden fermions is very 
small or the constituent fermion masses Mi and M 2 are accurately fine-tuned (or both), n° 
decays immediately into two ys. Since the stable DM particles can annihilate into two 7f°s 



Figure 9. DM annihilation into two ys. 
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Figure 10. DM annihilation into the SM particles via an internal S line and S — h mixing. In the actual 
calculation of the cross section we use the localized expression for the one-loop part, i.e. G v 2 g etc given in 
(37)-(39). 

with a huge DM conversion rate, there will be almost no DM left in the end. Since we want 
to assume neither a tinny Q nor accurately fine-tuned constituent fermion masses, we will 
not consider below the DM phenomenology based on the U(l)^, x U(l)§ flavor symmetry. 


2. (ii) SU{2) V x U(l)g 

Now we come to the case (ii) in (8), which means y\ — y 2 < yz- In this case the unstable 
NG boson is fj which can decay into two ys (and also into two S s if it is kinematically allowed). 
Because of SU(2) V , ft 0 is now stable and = m s o = m^± (see (21)). Furthermore, 
mfro and m^±, which are slightly larger than , are exactly degenerate in this case, i.e. 
m k ~ m k± = m k°■ H ie parameter region ( 20 ) we can further constrain the parameter 
space. Since y is a measure of the explicit chiral symmetry breaking and at the same time is 
the strength of the connection to the SM side, the smaller is y, the smaller is the DM mass, 
and the larger is the cutoff A. We have also found that for a given set of A h, A hs and A s 
the value (S) remains approximately constant as y varies, implying that ms also remains 
approximately constant because the Higgs mass ~ 126 GeV and ty = (h) — 246 GeV 
have to have a fixed value whatever y is. Consequently, the DM masses are smaller than 
ms, unless yi ^ 0.015 or A s and A hs are very small (or both). We have found, as long as 
we assume (20), that A s ~ 0.03 and A hs ~ 0.04 have to be satisfied to realize that S is 
lighter than DM. However, these values of A s and Xhs are too small for the DM annihilation 
cross sections into two Ss (diagrams in Fig. 8 ) to make the DM relic abundance realistic. 
In summary, there are three groups of DM in the SU(2)y case ( 8 ); the heaviest decaying 
SU(2) V singlet fj, two SU(2) V doublet ({K°,K~},{K + ,K 0 }) and lightest SU{2) V triplet 

IT 0 ) . 

Before we compute the DM relic abundance, let us simplify the DM notion: 

Xi = fj, X'2 to represent K°,K ± ,K°, y 3 to represent ft ± ,ft° (50) 

with the masses 

m\ = mfj, m 2 = m^ and m 3 = rri^ (mi > m 2 > m 3 ) , (51) 

respectively, where Xi are real scalar fields. There are three types of annihilation processes 
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which enter into the Boltzmann equation: 

Xi Xi ** X X , (52) 

Xi Xj ++ Xk Xi , (53) 

in addition to the decay of xi into two ys, where X stands for the SM particles, and 
the second process (53) is called DM conversion. There are two types of diagrams for 
the annihilation into the SM particles, Fig. 9 and Fig. 10. The diagrams in Fig. 10 are 
examples, in which a one-loop diagram and a tree-diagram are connected by an internal S 
or a S — h mixing. The same process can be realized by using the right diagram in Fig. 7 for 
the one-loop part. It turns out that there is an accidental cancellation between these two 
diagrams so that the velocity-averaged annihilation cross section is at most 10 -11 GeV -2 , 
unless near the resonance in the s- channel [51]. The effective 0-0-y interaction (23) can 
also contribute to the s-channel annihilation into the SM particles. However, as we have 
mentioned, the effective coupling Gk+k - 7 is very small in the realistic parameter space. 
For instance, m?^G K + K -^ f G K 25 ~ 10 —5 , where G K + K - 1 and G K 25 are given in (24) and 
(38), respectively. Note also that the DM conversion with three different DMs involved is 
forbidden by SU( 3)y. In the SU( 2)y x U(l)§ case, for instance, fj tt~ —> K° K~ is indeed 
allowed. However, it is strongly suppressed (G^n/G^^ ~ 1 CT 4 ), because SU(3)y is only 
weakly broken in the realistic parameter space. So we will ignore this type of processes, too, 
in the Boltzmann equation. 

Using the notion for thermally averaged cross sections and decay width (of xi) 

<vcr(ii;XX)> , <vcr(ir,jj) > , <T(1; 77 ) > , (54) 

the reduced mass 1 /p = (JTm” 1 ) and the inverse temperature x = g/T , we find for the 
number per comoving volume Y t = rii/s [ 66 ] 


dY 1 

dx 


= -0.264 g 1 / 2 


/iAfpL 


or 


{<va( 11-XX)> (F.F, -UU) 


+ <w(ll;22)>^riHi - 
xM Ph 


u 2 u 2 . 


Hi Hi + <w(ll;33)> - 


-0.602 g; 1 ' 2 
dY 2 


dx 


L ^ 

= -0.264 g 1 / 2 


Y 2 Y 2 

<F(1; 77) > (FI — Yi) , 

pMpL 


Y,Y, 


3 1 3 -r. 


Y 3 Y 3 


Y\Y\ 


(55) 


x z 


{<va( 22; XX) > (Y 2 Y 2 - Y 2 Y 2 ) 


Y 2 Y 2 , 


- <ucx(22; 11) > UU - t^xYY, 


Y 2 Y 2 


+ <w(22;33)>( Y 2 Y 2 - ^Y 2 Y 2 

13X3 


(56) 


dx 


= -0.264 g 


1/2 


pMpL 


x z 


{<u<t(33; XX) > (Y 3 Y 3 -Y 3 Y 3 ) 


y 3 y 3 . 


y 3 y 3 


- <va( 33; 11) > Y X Y X - - < w(33; 22) > Y 2 Y 2 - ~K 2 U 5 


Us > 3 , 


Y 3 Y 3 


(57) 


where g* = 115.75 is the total number of effective degrees of freedom, s is the entropy density, 
Mpl is the Planck mass, and Y % = rii/s in thermal equilibrium. Although XX)) is 

much smaller than 10 ~ 9 GeV -2 , we can obtain a realistic value of Qh 2 . The mechanism is the 
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following [ ]. If 2/3 does not differ very much from y\ = y 2 , the differences among mi, m 2 and 

m 3 are small. Then at finite temperature inverse DM conversions (which are kinematically 
forbidden at zero temperature) can become operative, because the DM conversions cross 
sections are large, i.e. 10 ~ 5 GeV 2 x phase space, as we have mentioned above. That is, the 
inverse conversion XsXz, X2X2 X1X1 ~t 7777 can play a significant role. 

The relic abundance Vth 2 is given by 


7 2 

\lih = 


Pc/h 2 


(58) 


where Y ioo is the asymptotic value of Y t , s 0 = 2890/cm 3 is the entropy density at present, 
p c = 3H 2 /&tiG = 1.05 x 10~ 5 h 2 GeV/cm 3 is the critical density, and h is the dimensionless 
Hubble parameter. Before we scan the parameter space, we consider a representative point 
in the four dimensional parameter space with Q — 1/3: 

y 3 = 0.00424, yi = y 2 = 0.00296, A s = 0.13, X HS = 0.06, X H = 0.135 , (59) 


which gives 

Qh 2 = (fti + 4D 2 + 3 Q 3 )h 2 = 0.119 , m s = 324.1 GeV, 

m 3 = mfj = 202.0 GeV, m 2 = m^ = 196.3 GeV, m 3 = m^ = 178.1 GeV, 

<va{ 11, 22, 33; XX) >= (9.29, 9.38, 1.26) x 10" 11 GeV -2 , 

< w(ll; 22) >= 4 < va(22; 11) >= 3.90 x 10“ 5 GeV“ 2 , 

< w(ll; 33) >= 3 < w(33; 11) >= 4.30 x 10“ 5 GeV“ 2 , 

<w(22;33)>= (3/4) <w(33; 22) >= 4.06 x 10" 5 GeV 2 , 


< T(l; 77 ) >= 6.45 x 10 " 13 GeV 1 , 

<na(ll; 77 ) >= 6.59 x 10 ' 14 GeV 2 = 7.73 x 10 " 31 cm 3 s ' 1 . 

Fig. 11 (left) shows Vth 2 (red), Qih 2 = Vt^h 2 (black), 4D 2 h 2 = 4flV 2 (green) and 3Vt 3 h 2 = 
3f \h 2 (blue) for the parameter values (59) as a function of the inverse temperature x = p/T. 
In Fig. 11 (right) we show the total relic DM abundance Q/ 1 2 as a function of yi(= y 2 ), where 
the other parameters are fixed as (59). Since a realistic value of Dh 2 for the SU(3)v case 
(9) can be obtained only near the resonance, i.e. ms/mDM — 2, the parameter space for the 
SU{2)v case ( 8 ) is considerably larger than that for the SU(3)v case (9). Note, however, 
that the realistic parameter space for the SU{2) V case is not continuously connected to that 
for the SU{3)v case, as we can see from Fig. 11 (right) (SU( 3)y means the point at yi = y 3 ). 


D. Indirect and Direct Detection of Dark Matter 

1. Monochromatic 7 -ray line from DM annihilation 

As we can see from Fig. 9, two DM particles can annihilate into two 7 s. Therefore, the 
charge Q of the hidden fermions can be constrained from the 7 -ray observations [68-70]. 
Since in the SU(2)y case the relic abundance of the 7 r dark matter is dominant, we consider 
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Figure 11. Left: The relic DM abundances, flh 2 (red), flih 2 = D^/i 2 (black), 4fl 2 h 2 = 4 fl^h 2 (green) 
and 3il 3 h 2 = 3 fl^h 2 (blue), as a function of the inverse temperature x = g/T for the parameter values (59). 
Though fj is almost in thermal equilibrium, its presence is essential for the K and tt numbers to decrease 
as x decreases. In the end the relic abundance of n dominates. Right: The total relic abundance flh 2 as a 
function of y\ (= y 2 ) for y 3 = 0.00424, where the other parameters are fixed as (59). 


here only its annihilation into two 7 s. We will take into account only the s-wave contribution 
to the annihilation cross section, and correspondingly we assume that p — p' — (m^, 0) 
and that the photon momenta take the form k = (m^,k) and k' = (m^, — k) with their 
polarization tensors e(k) = (0,e(k)) and e(k') = (0,e(k')) satisfying 

0 = e(k) ■ k = e(k) ■ k! = e(k) ■ p = e(k) ■ p 1 , 

0 = e(k') ■ k = e(k!) ■ k' = e(k') ■ p = e{k') ■ p , (60) 

respectively. 

To compute the annihilation rate we use the effective interaction (44). We find that the 
annihilation amplitude can be written as 

a b 

1 77 , ( 61 ) 

-tw 7 Z y J 

t 2 w zz 

where G n 2^2 is given in (46). Then (the s-wave part of) the corresponding velocity-averaged 
annihilation cross sections are 

a b 


r^tt/(u b'j — G n 2ry2 (k • k g^ h 


kjjK 


)x 


(va(7T7T —>■ a b)) = 


G^ml 


47T 


x < 


( 1 / 2 ) 


7 7 
7 Z 

(3/4)4(l-m|/ml) 1 /2 ^ 


4(l-m|/4ml) 


(62) 


The energy E y of 7 -ray line produced by the annihilation into 7 Z is mj( 1 — m|/4?n|). In 
practice, however, due to finite detector energy resolution this line cannot be distinguished 
from the E. f = rn^ line. Therefore, we simply add both cross sections. So we compute 
(vo) 11+1 z = (va(7T7r —> 77 )) + (va(nn —>■ 7 Z)) with Q = 1/3 as a function of for 
different values of A#, A# and Xhs , which is shown in Fig. 12 (right), where Vtti 2 is required 
to be consistent with the PLANCK experiment at 4a level [ ]. As we see from Fig. 12 
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Figure 12. Left: The Fermi Lat [69] (black) and HESS[ 68 ] (red) upper bounds on the velocity-averaged 
DM annihilation cross section for monochromatic 7 -ray lines, where this graph is taken from [ 68 ]. Right: 
The velocity-averaged DM annihilation cross section (vcr) 11 + 1 z as a function of m* with Q = 1/3. Since 
(t)<r) 77 + 7 z is proportional to Q 4 , our calculations can be simply extended to the case of an arbitrary Q. 
The red points are those for the SU{ 3)y case (9). 


(right) the velocity- averaged annihilation cross section is mostly less than 10 -29 cm 3 /s 
in the parameter space we are considering, and consequently the Fermi LAT and HESS 
constraints given in Fig. 12 (left) are well satisfied. The red points are those for the SU(3)v 
(9) case. 


The differential 7 -ray flux is given by 

, , dN*n 


dE~ 


oc (va)- 


(va 


)iz 


dN lz 

dE~ 


— {vir)- n+1 z <5 (£ 7 - inr, m) . 


(63) 


-^ 7 dEjry Ci-C / 7 Z 

Prospects observing such line spectrum is discussed in detail in [72-74], Obviously, with 
an increasing energy resolution the chance for the observation increases. Observations of 
monochromatic 7 -ray lines of energies of 0(100) GeV not only fix the charge of the hidden 
sector fermions, but also yields a first experimental hint on the hidden sector. 


2. Direct Detection of Dark Matter 


As we can see from Fig. 11 (left), the relic abundance of the K dark matter is about 
three orders of magnitude smaller than that of the 7 1 dark matter. Therefore, we consider 
only the spin-independent elastic cross section osi of 7 r off the nucleon. The subprocess is 
the left diagram in Fig. 13 (left), where • is the localized one-loop contribution (39), and we 
ignore the right diagram. The result of [ 3] can be used to find 

Z\ „ 0 fm N sin 29 ( 1 1 j m^rin 


°si = 


7T 


2v h mi 2 


mi 


mi 


(64) 


my + m-s- 

where G^s is given in (39), is the nucleon mass, and / ~ 0.3 stems from the nucleonic 
matrix element [ ]. We assume | cos 6 1 1 > 0.9 to satisfy the LHC constraint, where 9 is 
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m~ [GeV] 

Figure 13. Left: Subdiagrams contributing to the spin-independent elastic cross section asi off the nucleon. 
Since the relic abundance of K is negligibly small, the right diagram does not contribute. Right: The spin- 
independent elastic cross section a si of 7r as a function of m#. The red points are those for the SU(3)y 
case (9). The result should be compared with the XENON1T sensitivity of 0( 10 -47 ) cm 2 [77]. 


the h — S mixing angle. In Fig. 13 (right) we show in the m^-asi plane the area in which 
Vth 2 = 0.12 ± 0.01 (4cr) [ ] is satisfied. The predicted values of asi for ^ 150 GeV is 
too small even for the future direct DM detection experiment such as XENON1T, whose 
sensitivity is of O(10 -47 ) cm 2 [ ]. The smallness of asi results from the smallness of the 

coupling G n 2 S , whose smallness comes from small Yukawa coupling y\ and the accidental 
cancellation between the left and right diagrams in Fig. 7. The red points are those for the 
SU(3)v (9) case. We recall that the realistic parameter space for the SU(2 )y case is not 
continuously connected to that for the SU(3)v case, as one could see from Fig. 11 (right), 
in which y 1 = y 3 has to be satisfied for the SU(3)y case. 

If the relic abundance of the K dark matter were of 0(0.1), the non-zero K°-K°- r y/Z and 
K + -K~-'y/Z couplings shown in Fig. 13 would lead to a serious problem. Fortunately, the 
effective coupling is very small as we have already noticed: rri^Gx+K- 7 ~ 10 -6 , where this 
coupling for tt vanishes in the SU{2)y case. 

Note that because an accidental U{l)v (the hidden baryon number), not only the DM 
candidates, but also the lightest hidden baryons are stable. The hidden mesons in our model 
are neutral, while the charge of the hidden baryons b formed by three hidden fermions is 
Q b = 3 x Q. Let us roughly estimate the amount of relic stable hidden baryons and anti- 
baryons in the Universe, where we assume that the hidden proton and neutron are the 
lightest baryons in the SU(2) V case (8). As the hidden sector is described by a scaled- 
up QCD, the hidden meson-baryon coupling G^bb is approximately the same as in QCD, 
i.e. G^bb 13, which is independent of Q. Using this fact, we can estimate Y~ b = n b /s 
and obtain Y~ b ~ (0.4, 6, 9) x 10~ 16 for m b — 1, 5 and 8 TeV, respectively. There are 
severe constraints on Y~ b . The most severe constraints exist for Q b = 1, which come from 
the search of heavy isotopes in sea water [ ] and also from its influence on the large scale 
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structure formation of the Universe [ ]. We therefore conclude that Q b = 1, i.e. Q = 1/3 
is ruled out. Another severe cosmological constraint is due to catalyzed BBN [79], which 
gives Y~ b < 2 x 1CT 15 [80] (see also [81] ). The CMB constraint based on the Planck data is 
Q b h 2 ^ 0.001 [82], which can be satisfied in our model if m b ^ 6 TeV. (See also [83] in which 
the constraints in the Q^-DM mass plane are given, where these constraints are satisfied in 
a wide area of the parameter space of the present model.). 

In most of our analyses on DM here we have used Q = 1/3. The relic abundances of 
DMs depend on Q, because the decay rate of the neutral would-be DM y depends on Q. 
The change of Q can be compensated by varying the ratio of y 3 to y\ — 7 / 2 , as far as the 
difference of two hypercharges are not very much different. As for the indirect detection 
of DM, the annihilation cross section into two 7 s (62) being proportional to Q 4 should be 
multiplied with (3Q ) 4 for Q different from 1/3. The spin-independent elastic cross section 
(jsi (64) is independent on Q. This means that our basic results obtained in this paper can 
be simply extended to the case with Q different from 1/3. 


V. CONCLUSION 

We have considered a QCD-like hidden sector model [48-51], in which dynamical chiral 
symmetry breaking generates a mass scale. This generated scale is transmitted to the 
SM sector via a real SM singlet scalar S to trigger spontaneous breaking of EW gauge 
symmetry [48, 49]. Because the SM is extended in a classically scale invariant way, ’’Mass 
without mass” [1, 84] is realized in this model. Since chiral symmetry is dynamically broken, 
there exist NG bosons, which are massive because the coupling of S with the hidden sector 
fermions breaks explicitly the SU{rif)L x SU(nf)n chiral symmetry down to one of its 
diagonal subgroups. The mass scale of these NG bosons is calculable once the strength of 
this coupling and the scale of the QCD-like hidden sector are given. The smallest subgroup 
is the Cartan subalgebra U(l) nf ~ 1 . Because of this (accidentally) unbroken subgroup, the 
NG bosons charged under f/(l ) n/_1 are stable: There exist at least n 2 — nf DM candidates. 
We have restricted ourselves to n c = nf = 3, because in this case we can relate using 
hadrons the independent parameters of the NJL model, which we have used as a low-energy 
effective theory for the hidden sector. There are three possibilities: (i) f/( 1 )^, x U( 1 )^, (ii) 
SU{2)v x U(l)g, and (iii) SU(3)v, where the possibility (iii) has been studied in [51, 52], 
It turns out that the first case (i) is unrealistic, unless this case is very close to (ii) or (iii), 
or/and the hypercharge Q of the hidden fermions is tiny. This is because the lightest NG 
boson is neutral under U( 1)^, x so that it can decay into two 7 s and the stable DM 

candidates annihilate into them immediately. Therefore, we have mainly studied the case 
(ii) with yi = y-2 < V3- In this case the unstable NG boson is fj (the heaviest among the 
pseudo NG bosons) and can decay into two 7 s. The annihilation cross section into the SM 
particles via the singlet S is very much suppressed, except in the resonance region in the 
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s-channel annihilation diagram of DM. However, we have found another mechanism for the 
stable DMs to annihilate: If y 3 does not differ very much from y\ = y 2 , the differences 
among and rrifj are small. At finite temperature the inverse DM conversions (which 

are kinematically forbidden at zero temperature) can become operative, because the DM 
conversions cross sections are large ~ 10 -5 GeV~ 2 . Consequently, the realistic parameter 
space of the case (ii) is significantly larger than that of the case (iii), which has been obtained 
in [51, 52], 

With a non-zero Q the hidden sector is doubly connected with the SM sector; we have 
a bright hidden sector at hand. The connection via photon and Z opens possibilities to 
probe the hidden sector at collider experiments such as e + e~ collision [ ]. I 11 particular, 

the would-be DM, fj, can decay into two 7 s, which would give a smoking-gun event. 
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Appendix A: The NJL Lagrangian in the self-consistent mean field (SCMF) 

approximation 


Here we consider the NJL Lagrangian £njl (13) in the SCMF approximation [ ]. In the 

SCMF approximation one splits up the NJL Lagrangian (13) into the sum 

£njl = £0 + £-1 > (Al) 

where £/ is normal ordered (i.e. ( 0 |£/| 0 ) = 0 ), and £0 contains at most fermion bilinears 
which are not normal ordered. We find that £0 can be written as 

£0 = £/c + £d + £m , (A2) 

where 
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Here 77 ° stands for 0 O , and the meson fields are defined in (4). 


(AS) 


Appendix B: Determination of the NJL parameters G, Go and A 

As in [ ] we apply the NJL Lagrangian (A2) with g' = 0 to describe the real hadrons, 
where we assume SU(2)v and replace yiS by the current quark masses, i.e. yiS = y 2 S —)• 
m u , y 3 S —> m s . Then we compute the real meson masses m. n -,mK,m v ,m rt i and decay 
constants f n , Jk■ 

We obtain the following inverse meson propagators: 
r^±o 2 ) = r T o(p 2 ) 

= ~^G + + ( X _ 2ricI £ 2 ( p2 ’ m i’ m i) + §r n ^(m 3 ) , (Bl) 

rv±(p 2 ) = r A - 0 A 0 (p 2 ) 

= ~2 l ~G + + i 1 - ^§ ai ) ^(p 2 ,™!,^) + , (B2) 
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Parameter 

( 2G Q cd )- 1 / 2 

( g Q cd ) 1 / 5 

^QCD 

m u 

m s 

Value (MeV) 

361 

406 

930 

5.95 

163 


Table II. Values of the QCD NJL parameters obtained by fitting the pion and kaon decay constants 
and the meson masses, where we have assumed the SU(2)y flavor symmetry. 
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where the integrals /^ 2 (p 2 , m a , mf) and lf 2 {rn) are defined in appendix (C2), and 

mi = + <Ti - ^^(Ticr 3 , m 3 = m s + a 3 - ^^<7i 2 • 

The pion and kaon masses are the zeros of the inverse propagators, i.e. 

I\±(p 2 = mf 2 ) = 0 , T K ±{p 2 = m K 2 ) = 0 , 


(B6) 

(B7) 


while the 77 and rf meson masses are obtained from the zero eigenvalues of the real part of 
the r] 8 — rf mixing matrix. The wave function renormalization constants can be obtained 
from 


1 _ dV 7T ±(p 2 

Zj rrr 


Zk~ X = 


dV K ± 0 p 2 


dp 2 p2=m . 2 ’ dp 2 

and the pion and kaon decay constants are defined as 


< 0|TrV>7„7 5 l(<7i + s<r 2 )i/l|ir + (p) > = is/2f*p r , 

< 0|TrVi7„7 5 l((74 + ia^iplK*(p) > = isj2j K p v . 


(B8) 


(B9) 

(BIO) 


We use m,r, tur, m v , my, and Jr to determine the QCD NJL parameters. The best fit 
values of the parameters are given in Table II. I11 Table III we compare the meson masses 
and decay constants calculated in the NJL theory with the experimental values. As we see 
from Table III, the NJL r/ mass is about 16% smaller than the experimental value. This 
seems to be a general feature of the NJL theory [ ] (see also [86]). 
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Theory(MeV) 

Experimental value(MeV) 

UItt 

136 

140(tt ± ) 135(tt°) 

m-A' 

499 

494(A ± ) 498 {K°,K°) 

mn 

460 

548 


960 

958 

fn 

93 

92(tt-) 

Ik 

105 

HO(A-) 


Table III. Comparison of the NJL values with the corresponding experimental values. 


Appendix C: One-loop Integrals 


•Vacuum energy 

To compute the effective potential (18) we need the vacuum energy 

r d A k 

In det(f — 77i) 


him) = 


i( 2 tt) 


1 


A 4 In I 1 + ^~ 


A 2 \ 


m 4 In ( 1 + 1 + m 2 A 2 

rri 2 J 


(Cl) 


167T 2 y y A 2 
•Inverse propagator of dark matter 

There are two types of diagrams which contribute to the inverse propagator of dark matter: 


/^(p 2 ,m a ,m 6 ) = 


1$ (m) = 


d 4 l Tr(/ -^ + m a ) 7 5 (/ + m 6 )75 
i(2vr) 4 ((l -p) 2 -ml)(P - m 2 ) 
d A k TJi 1 


-7m 


A — TTi In ( 1-|-- 

m 2 J 


A 2 \ 


(C2) 


i(27r) 4 ( k 2 — 771 2 ) 167T 2 

These expressions are used to find DM masses and wave function renormalization constants 
given in (B8), respectively. 

• 0 - 0-7 amplitude 


1% Jp,p',m a ,m b ) 


= (- 1 ) 


dH TV(/ + m a ) 7 5 (/ - / + m b )Y(l + f + m h )75 


+ (p p',m a yy m 6 ), 


z(27t) 4 ((/ + p) 2 — m 2 )(/ 2 — m 2 )((/ — p') 2 — m 2 ) 

= -(pV y - p ,fl p u )(p + p')J<h 2 ( m ai m b )-\ - (C3) 

with p 2 = p' 2 , where • • • stands for higher order terms in the expansion of the external 
momenta, and 

1 1 
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The effective 0-0 -7 interaction Lagrangian is given in (23). 

• 0 - 7-7 amplitude 

The (j)(p)-'y(k)-'y(k / ) three-point function is needed to compute the decay fj into two 7 s 
(Fig. 4): 

I^(k, k', m) 

d 4 l Tr (/ — ft' + m) 7 M (/ + m)Y(t + ft + m) 75 


= (-l) 


i(27r) 4 ((/ + k ) 2 

1 7,/ _i_ ... 

c rvfy^rvo ^ 


m 2 ){l 2 — m 2 )({l — k') 2 — m 2 ) 


+ (fc O k 1 , p i/) 


(C5) 


47r 2 m “ 13 

The amplitude is thanks to 75 gauge invariant even for a finite A, i.e. 

k u I^ 2 (k, k', m) = k'^I^ik, k', m) = 0. The amplitude without 75 correspond to the 
S(p)- n f(k)-'y(k') three-point function, which we denote by I(f u s (k, k', m). This amplitude 
is not gauge invariant so that we need to apply least subtraction method [52] . The subscript 
0 indicates that the amplitude is unsubtracted, and we denote the subtracted gauge-invariant 
one by 7(0 (k, k'pm). In Appendix C we demonstrate how to use least subtraction method 
for this case. 


• 0 - 0 - 0-0 amplitude 

The 0(p)-0(p / )-0'(fc)-0'(A; / )- four-point function is needed to compute the DM conversion 
cross section (diagrams of Fig. 5 ): 
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At the lowest order in the expansion in the external momenta we obtain 
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(C12) 


where • • • stands for terms of 0(A 2 ) and higher. These expressions are used for the effective 
couplings defined in (32)-(35). 

•(j)-(j)-S amplitude 

To obtain the cj)(p)-(j)(p , )-S(k ) three-point function (Fig. 7) we need 
I^ s {p,p\rn a ,m b ) 

d 4 l Tr(/ + p + m 6 ) 7 5 (/ + m 0 ) 7 5 (/ - / + m b ) 
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At the lowest order in the expansion in the external momenta we obtain 
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= S ( _1 + ln ( A2 / m “) + '") ’ ( C16 ) 

/^2 S (0,0,m a ) = I$2 S (m a ) 

= 4 ^ ( A2 + 2m “ ~~ 3m “ ln ( A 7 m a) + ' " ) • ( C17 ) 

These expressions are used for the effective couplings defined in (37)-(39). 

•(jycjyS-S amplitude 
Similarly, 
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At the lowest order in the expansion in the external momenta we obtain 
/^ 2 S 2 ( 0 , 0 , 0 , 0 ,m a ,m 6 ) = 2I^ s2 (m a ,m b ) 
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— ln(A 2 /m^) + ■ • • ] , 
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771 

/^ 52 (0, m) = 2I^2 S 2 (m) = — (5 - 31n(A 2 /m 2 ) + • • •) . 

These expressions are used for the effective couplings defined in (41)-(43). 
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• 0 -<^>- 7-7 amplitude 

The next example is the <f>(p)-(f)(p')-'y(k)- r y(k') four-point function. The diagrams at the 
one-loop level are shown in Fig. 9: 

p\ k, k', m a , m b , m c ) 

_ [ d 4 l Tr(/ + m a ) 7 5 (/-/ + m b )Y{t +f ~ # + m c )Y(l + f + m b )^ b 

J i( 27t) 4 (Z 2 — m 2 )((Z — p') 2 — m 2 )((Z + p — k) 2 — 777 2 )((Z + p) 2 — ml) 

+ (p <H- p',k kp z/) T (A; A; 7 , p -H- 1 /) T (p p') , (C26) 
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= (-l) 


f(27r) 4 ((Z — fc ) 2 — m 2 )(Z 2 — m 2 )((Z _ p /^2 _ ml)({l + p — k) 2 — ml) 
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(C27) 


_ Tr(/ - #' + m) 7 ^(/ + m) 7 ^(/ + ^ + m) j 

j z(27t) 4 ((Z + k ) 2 — m 2 )(P — m 2 )((Z — Zd ) 2 — m 2 ) ’ 

The subscript 0 indicates that the amplitudes are unsubtracted, and therefore they are not 
gauge invariant. We apply least subtraction method to obtain gauge invariant amplitudes 
^, l!p 2 fiu an d , respectively. Since the realistic parameter space is close to that of the 

SU(2)y case ( 8 ), we consider them only in this case. At the lowest order in the expansion 
in the external momenta we obtain 

I^rik, k', m) + I*r(k, k\ 777) = (k ■ Per - W") + ■■■ (C 29 ) 
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& u (k,k',m) = 


1 


67r 2 m 


(k ■ k'g^ - Wk! v ) + 


(C30) 


in the large A limit. The result is used for the effective Lagrangian (44) and (61). 


Appendix D: Least Subtraction Procedure 


The cutoff A breaks gauge invariance explicitly and to restore gauge invariance we have to 
subtract non-gauge invariant terms from the original amplitude. In renormalizable theories 
there is no problem to define a finite renormalized gauge invariant amplitude. In the limit 
of A — y 00 the gauge non-invariant terms are a finite number of local terms, which can be 
cancelled by the corresponding local counter terms so that the subtracted amplitude is, up 
to its normalization, independent of the regularization scheme. To achieve such a uniqueness 
in cutoff theories, one needs an additional prescription. 

In [52] a novel method called “ least subtraction procedure ” has been proposed. The basic 
idea is to keep the subtraction terms to the minimum necessary. Consider an unsubtracted 
amplitude 

(A;p 1 ...pn s ,k 1 ...k ng ), (Dl) 

with n g photons and n s scalars (scalars and axial scalars). Expand the amplitude in the 
external momenta k's and p's: 

= E ■ < D2 ) 

m =0 

where Aq'^ /ln consists of m-th order monomials of the external momenta. I 11 general, 
Jlq 0 j 11 ,,, flns = (A; 0, • • • , 0) is non-vanishing and we can subtract it because it is not 

gauge invariant. We keep the tensor structure of as the tensor structure of the 

counter terms for A[™^ fln (m > 0) until a new tensor structure for the counter terms is 
required. We continue this until no more new tensor structure is needed. 


To illustrate the subtraction method we consider the S(p)-^(k)-^[k') three-point function, 


which is given by 

3 

•A 0 ^v(k,k') = ^ yjn c e 2 Q 2 

i —1 


dH Tr(/ - f + Mi)^ + MMt + f + Mi) 
i( 2vr) 4 ((l + k) 2 - M 2 )(l 2 - M 2 )((l - k') 2 - M 2 ) 


+(k tA k!, p o v) , 


(D3) 


where we use the on shell conditions k 2 = k' 2 = 0. Without lost of generality the amplitude 
can be written as 


A),nu(k, k') = A 0tg (k, k')g^ u + Ao tk (k, k')k^k' u + B 0tk (k, k')k v k'n . (D4) 

The last term does not contribute to the gauge invariance k u A^(k, k') = k'^A^k, k') = 0, 
and so we ignore it. The corresponding one-loop diagram is the one in Fig. 4 with fj replaced 
by S. According to least subtraction method, we expand the amplitude in the external 
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momenta k and k'. At the second order, for instance, we find 

4 2 }(k, k') = - n 4^(k ■ k') Y ,,. f ,f; A4 ,. f2 v, (2 A 2 + M? 


4n 2 
n c e 2 Q 2 


3M*(A 2 + M 2 ) 3 


i A 4 


(D5) 


VA 3M#nw (2A2+ nm 1 ' (D6) 

In the A —>■ oo limit the second order amplitude will be gauge invariant, but it is not at a 
finite A. Moreover, there are infinitely many ways of subtraction to make the second order 
amplitude gauge invariant. However, none of them is preferential. Least subtraction method 
uses the lower order amplitude, i.e. 

A 4 M* 




, k ') = 0 


(D7) 


(A 2 + M?) 2 ’ 

in this case, how to subtract the second order amplitude. At the lowest order in the derivative 
expansion, what is to be subtracted is unique; it is the g liu term. We keep this tensor 
structure as the tensor structure of the counter terms for higher order terms until a new 
tensor structure for the counter terms is required. However, in the case of Ao,nu(k, k') there 
will be no new tensor structure appearing in higher orders. This implies that Moy(/q k') 
remains unsubtracted (i.e. Ak(k,k') = Ao,k(k,k')) so that the subtracted gauge invariant 
amplitude is given by 

A 4 
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^ yi'n c e 2 Q 2 Mi 

-ZA- \9^k ■ k - k tl k v 
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dx 


where D = M 2 — 2xyk ■ k!. 
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pi — X 
'0 


dy 


2A 4 (1 - 4 xy) 

(A 2 + D 2 ) 2 D 2 


(D8) 
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